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We have studied the ground state phase diagram of a spin ladder with alternating rung exchange 
J" = Jx [1 + ( — in a magnetic filed, in the limit where the rung coupling is dominant. In this 
limit the model is mapped onto an XX Z Heisenberg chain in a uniform and staggered longitudinal 
magnetic fields, where the amplitude of the staggered field is ~ 5. We have shown that the magne- 
tization curve of the system exhibits a plateau at magnetization equal to the half of the saturation 
value. The width of a plateau scales as 5" , where = 4/5 in the case of ladder with isotropic 
antiferromagnetic legs and = 2 in the case of ladder with isotropic ferromagnetic legs. We have 
calculated four critical fields (-ff^i ™d corresponding to transitions between different magnetic 
phases of the system. We have shown that these transitions belong to the universality class of the 
commensurate- incommensurate transition. 

PACS numbers: 75.10.Jm Quantized spin models 



I. INTRODUCTION 



A theoretical understanding of the magnetic properties 
of quantum spin systems, in particular of spin 5 = 1/2 
isotropic antiferromagnetic two-leg ladders, has attracted 
a lot of interest for a number of reasons. On the one hand, 
there was remarkable progress in recent years in the fab- 
rication of such ladder compounds^ On the other hand, 
spin-ladder models pose interesting theoretical problems 
since antiferromagnetic two-leg ladder systems have a gap 
in the excitation spectrum and, in the presence of a mag- 
netic field, they reveal an extremely rich behavior, domi- 
nated by quantum effects. These quantum phase transi- 
tions were intensively investigated both theoretically'^^ 
and experimentally. 

The discovery of a magnetic field induced gap in the 
Cu-bensoate2i and other spin chain materials^ have in- 
creased the interest for magnetic quantum phase tran- 
sitions which are determined by the combined effects of 
the uniform and staggered components of the effective 

magnetic field,^3,24,25,26,27,2S,29,30 

In this paper we study the ground state magnetic phase 
diagram of the spin S* = 1/2 two- leg ladder with alter- 
nating rung exchange given by the Hamiltonian (Fig. ^ 

+ j±^[i + (-ir<5]s„,i-s„,2, (1) 

n 

where S„^q is a spin S' = 1/2 operator of rung n 
(n=l,...,N) and leg a [a — 1, 2). The interleg coupling is 
antiferromagnetic, = J_l(1 ± <5) > 0. 

We restrict our consideration to the limit of strong 
rung exchange 3> |J|||,(5J_l and map the model onto 
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FIG. I: The ladder with alternating rung exchange 



spin- 1/2 XXZ Heisenberg chain in the presence of both 
longitudinal uniform and staggered magnetic fields, with 
the amplitude of the staggered component of the mag- 
netic field proportional to '-^ SJ±. We study the ground 
state phase diagram of the effective spin-chain model and 
show, that the alternation of the rung-exchange leads to 
the dynamical generation of a new energy scale in the 
system and to the appearance of two additional quantum 
phase transitions in the magnetic ground state phase dia- 
gram. These transitions manifest themselves most clearly 
in the presence of a new magnetization plateau at mag- 
netization equal to one half of its saturation value (see 
FigEJ- The magnetic phase diagram is characterized by 
the following four critical fields: the field H~i, which cor- 
responds the the transition from a gapped rung-singlet 
phase to the gapless paramagnetic phase; the critical 
fields H^^ and H~2 which mark end-points of the magne- 
tization plateau and the saturation field i/^ . The width 
of the plateau scales as S'^ , where ly = 4/5 in the case of 
a ladder with isotropic antiferromagnetic legs and ly — 2 
in the case of a ladder with isotropic ferromagnetic legs. 
Therefore this magnetic phase diagram is generic for a 
standard isotropic ladder with alternating rung exchange. 
However, in the case of a ladder with ferromagnetic legs 
and frustrating diagonal interleg exchange, the interme- 
diate magnetization plateau dissappears for sufficiently 
strong ferromagnetic diagonal coupling. 
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FIG. 2: Schematic drawing of the magnetization (in units of 
saturated magnetization Maat) of a two- leg isotropic ladder 
with alternating rungs as a function of the external magnetic 
field. 



II. DERIVATION OF THE EFFECTIVE 
HAMILTONIAN 

In this section we derive the effective spin-chain 
model to describe the strong rung-exchange hmit J± 3> 
((5Jj^), I J|| I of the model (^). To obtain the spin chain 
Hamiltonian we follow the route already used to study 
the standard ladder models in the same limit of strong 
rung exchange!^ We start from the case Jy = 0. In 
this limit the system decouples into a set of noninter- 
acting rungs with couplings and JJ. In this case, 
an eigenstate of Ti. is written as a product of rung 
states. At each rung two spins S„^; and S„.2 are ei- 



ther in a singlet state |s° 



1 

V2 



(I Ti) 



it)) or in 

one of the triplet states: |t+) = 111), = | H) and 
|t°) = ^(1 Ti) + I iT))- Their energies are respectively: 
E{s°) - -3J^/4, E{t+) = Jl/A-H, E{tl) = and 
= J«/4 + H, where J_l[1 + (-1)"(5] 

When H is small, the ground state consists of a prod- 
uct of rung singlet. As the field H increases, the energy 
of the state |tj„_i) decreases and at H — Hd = JJ this 
state is degenerate with \s 



2n-l/ 



Thus, at iJ = Hci the 
ground state of a given odd rung undergoes a transition 
from the singlet |s2„_i) to the triplet |i^„_i) and the 
total magnetization of the system jumps discontinuously 
from zero to 0.5Msat = N/2. With further increase of 
the magnetic field, at Hd < H < Hc2 the magnetiza- 
tion remains constant. However since for Hc2 > the 
energy of the state |tj„) is lower then |s2„), the magneti- 
zation once again increases discontinuously from O.bMgat 
to Msat ^ N for H^2 = Jt- 

For J|| 7^ these abrupt transitions are broadened 
into intervals < H < H^^ and iJ^ < H < i?^, 
respectively. Two different scenarios are possible. Ei- 
ther H^^ < H~2 and the magnetization plateau with 
M = 0.5Msat remains. Or > i/^ and alternation 
of the rung exchange is irrelevant. In the latter case, the 



model shows a similar behavior as the standard two-leg 
ladder in a magnetic field: with increasing magnetic field 
in the range between and i?^ the magnetization con- 
tinuously evolves from a nonmagnetic phase at H < H^^ 
into the fully polarized ferromagnetic state at H > iJ^. 

The easiest way to obtain the effective model is to split 
the Hamiltonian into three parts: 



H = H, 



H, 



Hi, 



H, 



(o) 



H 



and 



Hi, 



- Jj_ "* S2„-l,l • 82,1-1,2 

n 

- Hcl 'S'Ira-l.a i 

n,a 
N/2 

- </j_^'' 82™, 1 • S2n,2 

n=l 

- iic2 '5'2ti,q 

- t^ll ^ ^ ^ qS„^i et 

n a 

- (H - Hcl) S2n-l,a 

- {H - Hc2) S2n,a ■ 



(2) 



(3) 



(4) 



The ground state of TCq is 2^ times degenerate, since each 
rung can be in the state | sq) or \t^) and the first excited 
state has an energy of the order of J± . Hint will lift the 
degeneracy in the ground state manifold, leading to an 
effective Hamiltonian that can be derived by standard 
perturbation theory.^ 

Let us start by introducing pseudo-spin r = 1/2 oper- 
ators, T„ which act on these states as 



tJ So >n 
rn\ So >r 



1, 

So >n , 
\t+>n, 

0, 



1, 



\t+ >n , 



r+|t+>„ = 0, (5) 

Tn\t+ >n= \so>n ■ 



The relation between the real spin operator S„ and the 
pseudo-spin operator r„ in this restricted subspace can 
be easily derived by inspection. 



(6) 



Using 0, to first order and up to constant, we easily 
obtain the effective Hamiltonian 

Heff = Y.^Jxy («+i + r^T^+i) + J.t:X+i} 

n 



3 



where 



J xy 


= ^11 , Jz 




(8) 




= H-Jj_- 


J|l 

" Y' 


(9) 




= SJj_. 




(10) 



A similar argument can be used to determine -ff^^. It 
is useful to make a particle-hole transformation and esti- 
mate from the condition where the transformed hole 
band starts to fill. This gives 

H+ =J^ + J„ + + {SJ^y . 



(16) 



Thus the effective Hamiltonian is nothing but the XX Z 
Heisenberg chain, with anisotropy Jz/Jxy = A = 1/2 in 
uniform and staggered longitudinal magnetic fields. It is 
worth to notice that a similar problem has been studied 
intensively in recent yearsi^^i^^'^^i^^'^'^'i^^'^^i^" 



III. MAGNETIC PHASE DIAGRAM 



The first critical field H^^ and the saturation 
field 



We first calculate of the critical field , correspond- 
ing to the transition from a gapped rung-singlet phase 
to a gapless paramagnetic phase, and the saturation field 

For H < H^^ the ground state of the system corre- 
sponds to the gapped rung singlet phase with zero mag- 
netization. For H > the system is in the fully po- 
larized ferromagnetic phase. The easiest way to express 
-ff^ and terms of ladder parameters J|| , J± and 

S is to perform the Jordan- Wigner transformation which 
maps the problem onto a system of interacting spinless 
fermions: 

^^f = i^(a+a„+i + /i.e.) + V^^p„p„+i 



^[Aio + (-l)>i]p„ 



where 



^=2^11 



Mo = 7;J\\ + hlff , Ml = 



(11) 
(12) 

(13) 



The lowest critical field corresponds to that value of 
the chemical potential /xqc for which the band of spinless 
fermions starts to fill up. In this limit we can neglect the 
interaction term in Eq. (|ll|l and obtain the model of free 
massive particles with spectrum 



(14) 



The chemical potential corresponding to H^^ is given by 



nl, i.e. 



(15) 



B. Magnetization plateau: H^^ and H^2 

To determine the values of the remaining two crit- 
ical fields and H~2 we consider the model Qfor 



For hOf^ 



the Hamiltonian iQ with 
anisotropy parameter |A| < 1 is known to be critical. 
The long wavelength excitations are described by the 
standard Gaussian theory with Hamiltonian^ 



K 



leg 



(17) 



Here (j){x) and 6{x) are dual bosonic fields, dt(j) — VsdxO, 
and satisfy the following commutational relation 



{x),9{x)]^t/2. 



(18) 



The velocity of spin excitation Vs is fixed from the Bethe 
ansatz solution as 



K 



■sin {■k/2K), 



2K - 1 

where the spin-stiffness parameter K is given by 
K 



2 1- 



arccos 



A) 



(19) 



(20) 



Thus the parameter K increases monotonically along the 
XX Z critical line — 1 < A < 1 from its minimal value 
iiT = 1/2 at A = 1 (isotropic antiferromagnetic chain), 
to unity at A = (the XY chain) and diverges at the 
ferromagnetic instability point of a single chain A = — 1 . 

To obtain the continuum version of the Hamiltonian 
O we use the standard bosonization expression of the 
spin operator— 

^ dx4> + (-1)" - sin(\/4^0) . 

TT TT 

where A is a non-universal real constant of the order of 
unit}^, and get the continuum Hamiltonian 

Hbos = I dx[^[{dx(j^f + {dx6f] 

■sm{VI^(l>) - hlffJ —dxA . (21) 

The Hamiltonian H21|) is the standard Hamiltonian for 
the commensurate-incommensurate transition which has 
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been intensively studied in the past using bosonization-^ 
and the Bethe ansatzJi Below we use these results to 
describe the magnetization plateau and the transitons 
from a gapped (plateau) to gapless paramagnetic phases. 

Let us first consider /i^j^ = 0. In this case the contin- 
uum theory of the initial ladder model in the magnetic 
field H = J_L-|- J||/2 is given by the quantum sine-Gordon 

(SG) model with a massive term ~ ^\ff sin(-\/ AnKcj)) . 
From the exact solution of the SG model2& it is known 
that the excitation spectrum is gapless ioi K > 2 and 
has a gap in the interval < K < 2. At K = 1 the 
sine-Gordon model is equivalent to the theory of free 
massive fermions with m = ft-gj^-- At 1 < X < 2 the 
excitation spectrum of the model consists of solitons and 
antisolitons with mass M, while for < if < 1 the 
spectrum contains also soliton-antisoliton bound states 
("breathers"). The exact relation between the soliton 
mass M and the bare mass h^jr^ = SJ± is given hy^ 



M = JfiK) {SJ±/J\\) 



l/(2-K) 



where 



CiK) = 



mi 



2v' 



r(i - K/2) 

2T{K/2) 



l/{2-K) 



(22) 



(23) 



It is straightforward to get from ()20(l . that at A = 1/2 
the spin stiffness parameter K = 3/4. Therefore, for 
h^eff ~ ^ sine-Gordon Hamiltonian (|21|l is in the 
strong coupling (massive) regime. In this case the low- 
energy behavior of the system is determined by the 
strongly relevant staggered magnetic field (i.e. alternat- 
ing part of the rung exchange), represented by the term 
sin(-\/37r(/)). In the ground state the field <j) is pinned 
in one of the minima of the staggered field potential 



(24) 



377010) = -7r/2 + 27rn. 



In view of H21|l we conclude that this state corresponds 
to a long-range-ordered antiferromagnetic phase of the 
effective Heisenberg chain ((TJ, i.e. to a phase of the initial 
ladder system, where odd rungs have a dominant triplet 
character and even rungs are predominantly singlets. 

At h^^ff 7^ (i.e. H ^ J± + J\\/2) the very presence of 
the gradient term in the Hamiltonian (|21|l makes it nec- 
essary to consider the ground state of the sine-Gordon 
model in sectors with nonzero topological charge. The ef- 
fective chemical potential ~ /i^j^ \J^dx4> tends to change 
the number of particles in the ground state i.e. to create 
a finite and uniform density solitons. It is clear that the 
gradient term in (|21|l can be eliminated by a gauge trans- 
formation (j) — > (ps+h^ff^J^x, however this immediately 
implies that the vacuum distribution of the fieled (j) will 
be shifted with respect of the minima (|24|l . This compe- 
tition between contributions of the smooth and staggered 
components of magnetic field is resolved as a continuous 
phase transition from a gapped state at \h^ff\ < M to a 



gapless (paramagnetic) phase at \h^eff \ > M, where M is 
the soliton massi^ 

For our effective Hamiltonian {Tj) with A = 1/2, the 
spin stiffness parameter K is 3/4 model fEa. (|20|l ') and 
the commensurate-incommensurate transition in the ef- 
fective sine- Gordon theory at ±^e// — ^ gives two ad- 
ditional critical values of the magnetic field 



Ji + J,l/2- J,|Co {6J±/J\\) 



4/5 



and 



Hc2 = J± + ^11 /2 + J||Co {5J±/J\\) 



4/5 



(25) 



(26) 



where Co = C(3/4) = 1.11428. 

As usual in the case of quantum commensurate- 
incommensurate transition transitions, the magnetic sus- 
ceptibility of the system shows a square-root divergence 
at the transition points: 



H) 



X{H) = 



h-h: 



-1/2 



-1/2 



for H < H+ 
for H+ <H< 

for H > HJ, 



C. Magnetic phase diagram 

Summarizing the results of the previous subsections we 
obtain the following magnetic phase diagram for a lad- 
der with alternating rung exchange (see Fig. [2J). For 
H < the system is in a rung-singlet phase with 

zero magnetization and vanishing magnetic susceptibil- 
ity. For H > H~i some of singlet rungs melt and 

1/2 

the magnetization increase as (if — ii^i) ■ With fur- 
ther increase of the magnetic field the system gradu- 
ally crosses to a regime with linearly increasing magne- 
tization. However, in the vicinity of the magnetization 
plateau, for H < H^^ this linear dependence changes and 
the magnetization once again shows a square-root behav- 
ior M- \Msat {H+^ - Hf'^. For fields in the inter- 
val between H'^i and H~2 the magnetization is constant 
M — O.hMsat- At H > H~2 the magnetization increases 

1 /2 

as M ^ Q.^Msat + (ii — H^2) ' then passes again 
through a linear regime until, in the vicinity of the satu- 



ration field ii^, it becomes M ^ M^at 



1/2 



The width of the magnetization plateau at M 
0.5Msat is given by 



h: 



H+ ^ 2J,| (<5J^/J,i) 



4/5 



(27) 



D. Ladder with ferromagnetic legs 



The existence of a magnetization plateau at M = 
0.5Msat is not limited to the ladder with antiferromag- 
netic exchange, but is also found for a ladder with 
isotropic ferromagnetic legs (Jy — — |J||| < 0) coupled 



5 



by an antiferromagnetic rung exchange {J± > 0). In this 
case the effective spin chain model is also given by the 
Hamiltonian 0, but with different parameters: 



kill, 
5Jj_ . 



(28) 

(29) 
(30) 



Thus, in the case of ferromagnetic legs the anisotropy 
parameter of the effective XX Z is A — —1/2 and con- 
sequently the spin stiffness parameter K is given by 
K — 3/2. The equivalent sine-Gordon theory (|21|l with 
a massive term ~ /i^j^ sin(-\/67r0) remains in the gapped 
strong-coupling regime. Using Eq. (|22|l we find an ex- 
citation gap M ~ |J|[| ((5Jj_/|Jj||) . Correspondingly the 
width of the magnetization plateau in this case equals 



- H+ 0, 2J|| (W^/J|| 



IV. GENERALIZED LADDER 



(31) 



In this section we consider a generalized ladder model 
with frustrating (diagonal) interleg interactions. The 
Hamiltonian is given by 

a. 

n,a n,a 

+ Jj_ ( Sn,l • ^n-\-l,2 + ^n,2 ' S^^+i j ) 

n 

+ J±5I(l + (~l)"'^)S".i-S",2. (32) 

n 

Lat us first consider the case of antiferromagnetic legs 
(J|| > 0). Assuming J± » J||, \J'j_\,5J± we easily obtain 
the parameters of Hef / , 



2 Jii - ' 



H -J_i_ 



J|| 

Y 



2 



(33) 

(34) 
(35) 



Below we briefly discuss following limiting cases: (i) K ~ 
1/2 (A 1), (ii) = 1 (A = 0) and (iii) K ^ 2 {A = 
-1/V2). 

The case (i) is reached for antiferromagnetic diagonal 
exchange J'j_ — 1/3 Jy. In this case the width of the 
magnetization plateau is 



- H+ 0, 2J|| (<5J^/ J|| 



\2/3 



(36) 



The case (ii) is reached for ferromagnetic diagonal ex- 
change = — ^||, where the effective continuum theory 



becomes the theory of free massive fermions. At K 
the spin gap at h^^j = is 

M = 5J±_ . 



(37) 



Correspondingly, the width of the magnetization plateau 
also scales linearly in 5. 

The frustrating diagonal exchange can change qualita- 
tively the magnetic phase diagram of our model only in 
the rather special case of a ladder with ferromagnetic legs 
and ferromagnetic diagonal exchange. In this case the 
anisotropy parameter of the effective spin-chain model is 
given by 



1 J|| - Jj 
" 2 J|| + J'^ 



(38) 



This relation implies that for 



A<J'l 



\/2-l 
V2 + 1' 



(39) 



the spin stiffness parameter K exceeds 2. In this pa- 
rameter range the sine-Gordon model is gapless and the 
magnetization plateau disappears. 



V. CONCLUSION 

We have studied the ground state phase diagram of 
a spin S = 1/2 two- leg ladder with alternating rung- 
exchange J" = Jj^ [1 + (— 1)"(5] in a magnetic field. We 
have shown that in a wide parameter range the magne- 
tization curve exhibits a plateau at one-half of its satu- 
ration value. The width of the plateau, is proportional 
to the excitation gap in the system at M = 0.5Msat and 
scales as S'^ . The critical exponent has a value = 4/5 in 
the case of a ladder with isotropic antiferromagnetic legs 
and !/ = 2 in the case of a ladder with isotropic ferromag- 
netic legs. We have also shown that in the ladder with 
frustrating diagonal interleg exchange the plateau effect 
is stronger and for realistic values of diagonal exchange 
the critical exponent reaches the value v = 2/3. 

We have also shown that in the case of ladder with fer- 
romagnetic legs and with strong diagonal ferromagnetic 
intraleg interaction, the magnetization plateau M — 
0.5Msat is absent. 

To conclude, we briefly comment the possible sponta- 
neous appearance of an alternating rung exchange as a 
spin-Peierls instability. Assuming that in the harmonic 
approximation the lattice deformation energy per rung is 
given by Edef ~ and estimating the magnetic conden- 
sation energy as Emagi^) — Emag{0) ~ —5^'^, we conclude 
that such an instability is possible for an antiferromag- 
netic ladder for a magnetization M = Q.5Msat, i-e. for 
H = Hc~ Jj_ -h0.5J||. 
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